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Abstract 



The negatively charged quasi- two-dimensional exciton {X~) is studied 
numerically in the presence of a uniform perpendicular B-field and various 
in-plane confinements. The corresponding photoluminescence (PL) spectra, 
including shakeup processes, are calculated. These calculations provide a 
quantitative explanation of recent experimental PL spectra. The photolumi- 
nescence signatures of charged excitons in quantum dots are also examined. 
Electron-hole pair distribution functions are calculated which show the struc- 
ture of the charged exciton states (both ground and excited states). Increasing 
excitation reflects a competition between two effects: partial ionization of ei- 
ther one or two electrons. 
PACS: 78.20.Ls 78.66.-w 73.20.Dx 
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I. INTRODUCTION 



Excitons (X) are crucial to the understanding of optical spectra from quantum wells with 
a low electron density. If the electron density is increased the excitonic effects are destroyed 
by phase-space filling from the other electrons and a Fermi Edge Singularity (FES) may 
be seen. Recent experiments have shown that as the electron density is lowered, the FES 
evolves into a discrete line which is not an exciton but is instead attributed to the negatively 
charged exciton {X~ = X + e^) [jT| [Q]. The negatively charged exciton was first predicted 
by Lampert in analogy with the various ions of hydrogen |Q. Further variational studies 
have been carried out in both three dimensions (3D) [^] and two dimensions (2D) [^]. It is 
found that the restriction to 2D increases the binding energy of the second electron to about 
ten times its 3D value. This should facilitate the observation of the X~ and, consequently, 
most experimental studies are carried out in the quasi-2D environment of the quantum well. 
Negatively charged excitons were first observed in CdTe/Cdi_xZnxTe quantum wells 
and were later seen in GaAs/Gai_xAlxAs § g H lTl|]. Most of the experimental 
spectra are obtained using polarized photoluminescence (PL) techniques. A constant B-field 
is often applied perpendicular to the well since this allows a charged exciton to form even at 
electron densities for which an FES would be seen in zero field |0 . Recently the low energy 
tails of the PL spectra have been examined and evidence of shakeup processes have been 
found ||12|. These occur when the charged exciton undergoes electron- hole recombination 
leaving the remaining electron in an excited state. Positively charged excitons, consisting 
of an electron and two holes {X~^ = X + h'^) have also been observed in p-tpye quantum 
wells [| [|ll |4|. 



Recent theoretical studies of the charged exciton complex consider ideal 2D systems 15 



and are often restricted to the lowest Landau level [|T6|. In a previous paper we showed that 
the combined effects of (a) higher Landau level mixing and (b) the form factor emulating 
the finite quantum well width (i.e. quasi-2D system), are essential for understanding recent 
experimental X~ spectra P 
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In this paper accurate numerical solutions are presented for a quasi-2D negatively charged 
exciton (X^) consisting of two electrons and a hole situated in an in-plane confinement 
potential and a perpendicular uniform B-field. The properties are studied over a wide 
range of experimentally accessible B-fields and in-plane confinements. We consider weak, 
intermediate and strong in-plane confinements, and use a range of initial angular momenta 
for the charged exciton. We calculate the development with B-field of the following: 1) The 
photoluminescence spectrum 2) the shakeup lines in the PL spectrum and 3) the electron- 
hole pair distribution functions. For the case of a weak in-plane confinement we can compare 
our results to recent experiments in high mobility quantum wells. We find that theory and 
experiment match very closely in both the main spectrum (Sec. |T|) and the shakeup lines 
(Sec. For stronger confinements our results should be compared to experiments carried 
out in quantum dots. To our knowledge no charged exciton species have yet been knowingly 
observed in quantum dots. We predict rich signatures characterizing the X~ in the quantum 
dot and hope that these results may stimulate further experimental study in such systems 
(Sec. fT| and Sec. |IV|). Finally our study of the electron-hole pair distribution functions 
allows us understand the structure of the charged exciton species and how this relates to 
the features of our calculated spectra (Sec. |V|). We also analyse the distribution functions 
of the excited X~ states finding that increasing excitation refiects a competition between 
two effects: partial ionization of either one or two electrons. 

II. MODEL 

The negatively charged exciton consists of two electrons and a hole interacting in a 
constant magnetic field according to the following Hamiltonian: 
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where rrie — 0.067mo and nih — 0.48mo are the masses of the conduction electrons and 
heavy holes respectively for GaAs. If the interactions Vee, and Veh are set to zero and the 
symmetric gauge is chosen (i.e. A(r) = (B x r)/2) the electron and hole single particle 
states are as follows: 
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is ^ normalisation, LI^I are associated Laguerre polynomials, tiuje and fiujh are the 
parabolic confinements for the electron and the hole respectively, and uj^ f^ — eB/nich is 
the electron/hole cyclotron frequency. Also uJe,h = ((f^e/i)^ + ^^Ih)- The single particle 
wavefunctions have characteristic length scales l^ f^ — h,/me,hWe,h- In this symmetric gauge 
the single particle wavefunctions have two quantum numbers, n and M. In the limit of weak 
confinement (i.e.(c<jg ^ -iiue^hY) the magnetic field dominates and the single particle 
states form Landau levels. For electrons, states with Me < and n=0 are degenerate 
forming the lowest Landau level; raising n or making Mg positive raises the electron energy 
and puts it in a higher Landau level. For the holes the positive states with n=0 form 
the lowest hole Landau level. 

In the case of strong confinement {i.e.{ijjl f^'^ <S 4(a;e,?i)^) the single particle states are 
essentially oscillator states of the parabolic confining potential, where aJe,/i '2LJe,h- For 
a given n, states with the same \M\ are hence degenerate. Raising either n or \M\ raises 
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the energy by an integral number of hue^h- For both electrons and holes, the states have 
definite angular momentum, hMe^h, arising from the rotational invariance of the system. 
The experimental data referred to in this paper is all for weak confinements. For simplicity, 
therefore, we refer to Landau levels throughout this paper noting that the Landau level is 
no-longer degenerate when the confinement is strong. The position of the particle in a given 
single particle state is closely related to the angular momentum of the state, M. For particles 
in the lowest Landau level (i.e. n=0): 

(r^) = 2f (|M| + 1). (5) 

This paper considers three confinements: weak {i.e.{u;lf^y ^ 4(a;e,/i)^), intermediate 
{i.e.{uj2fi)'^ ~ 4(a;e,/i)^) and strong (i.e.(a;g ^j)^ <S 4(a;e,/i)^). The cyclotron energies are 
hcu^ = 17.2meV and hcuf^ = 2.4meV at B=10 T. The confinements are chosen to be huje =1-75 
meV in the weak confinement limit, hue —5.58 meV in the intermediate confinement limit 
and hcUg =17.6 meV in the strongly confined limit. For ease of comparison we chose u)h such 
that LJe/'^h = 'tTT'h/^e, thereby ensuring that the electron and hole wavefunctions have equal 
characteristic lengths, l^ j^ — h/rrie^hCJe^h — l^- In the case of weak confinement the magnetic 
field dominates and we can compare our results to those obtained from high mobility quan- 
tum wells. For intermediate and strong confinements the charged exciton is effectively in a 
quantum dot. In the case of intermediate confinement, the B- field dominates at reasonably 
large B while the confinement dominates at low B. The crossover occurs when cu^^ '-^ '^i^e,h 
i.e. when B~6.5T. As expected we find interesting changes of behaviour around the crossover 
field strength. 

The interaction potentials are Coulombic. In order to obtain a tractable model it is 
assumed that the single particle wavefunctions separate: \E'e = 0e(^e)'?/'e(re) and \E'/i = 
<Ph{zh)'4^h{T^h)- A specific form for 0e,fe is chosen, thereby freezing out the z-motion and 
yielding a quasi-2D model. There are two different limits we can choose for 0e,/i- In one 
limit the electrons and holes occupy the same plane but have their charge distributions 
smeared in a rod- like distribution along the perpendicular direction (rod geometry). In 

5 



the other hmit the electrons and holes exist on separate planes separated by a distance d 
(biplanar geometry), their charge distributions being unsmeared. The biplanar geometry 
is the relevant limit of systems such as a quantum well in a polarizing electric field which 
separates the electrons and holes. The rod geometry is suitable for quantum wells in which 
the electron and hole distributions are not separated in the direction perpendicular to the 
well. We showed in an earlier publication that the biplanar geometry doesn't exhibit 
a PL signature of charged excitons for d> I. This lack of X~ signature has been confirmed 
experimentally for a quantum well in a perpendicular electric field. Consequently we 
concentrate on the rod geometry in this paper, for which: 

^eA^e,.)? = (6) 

otherwise 



In this work we use a rod length L = 31. At B=10 T, L = 240A for weak confinement, 
L = 220A for intermediate confinement and L = 170A for strong confinement. 

Before proceeding with any calculation of charged exciton states it is important to iden- 
tify the symmetries of the full Hamiltonian as these lead to well-defined quantum numbers, 
thereby allowing us to work within a restricted Hilbert space. For the charged exciton the 
Hamiltonian in the symmetric gauge is rotationally invariant and also conserves the total 
spin. The good quantum numbers for the charged exciton are therefore the total orbital 
angular momentum hM, the total spin hS{S + 1) and the total z-component of the spin, 
hSz- Spin conservation, together with the fact that the wavefunction must be antisymmetric 
under interchange of the two electrons, allows us to label the charged exciton wavefunctions 
as either triplet or singlet. The singlet states have S=0 and = with an orbital wave- 
function which is symmetric under electron interchange. The triplet states have S=l and 
Sz = ±1 or with an orbital wavefunction which is antisymmetric to electron interchange. 
Thus the charged exciton wavefunctions have the form: 

\E'(M, S = 0, Sz = 0) = ^ ^niMi,n2M2,nhMh^niMi,n2M2,nhMS^^ 

niAIi >n2M2 
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x{s = o,s, = o) 

'^{M,S = l,Sz = ±1,0) = ^ C'^iAfi,n2Af2,nhMh0niAfi,n2M2,nhMh(^) 

niMi>n2M2 

x{S=l,Sz = ±l,0) (7) 

where 

0niAfi,n2Af2,nhMh(^) = — {0niAfi (ri)0n2Af2 (''2) + 0ni A/i (r2)0n2M2 (fl) } 0n^Mh (fh) 



^niMi,n2M2,nhMh 



v/2 

(^) = {0niAfi(ri)</)„2Af2(r2) - 0niAfi(r2)0„2Af2(ri)}0n,„A4(rh) 



Here x('S' = 0, Sz = 0), x{S = l,Sz = ±1, 0) are the usual singlet and triplet spin states re- 
spectively. The total angular momentum M = M1 + M2 + M/J. The restriction riiMi > n2M2 
is required to prevent overcounting of states. In the particular limit that the basis is re- 
stricted to the lowest Landau level, the interaction Veh = Vee and the in-plane confinement 
is zero, there is a well documented hidden symmetry This forces the photolumi- 



nescence energy from a triplet charged exciton to be indentical to that of the lowest exciton 
line. In our system this symmetry is broken by the in-plane confinement and by Landau 
level mixing, the latter being crucial in the case of weak confinement ||T^. This breaking of 



the hidden symmetry is vital if any signature of the triplet charged exciton is to be seen in 
photoluminescence . 

The effects of the B-field and the Coulomb interaction compete: from a perturbation 
theory perspective the expansion parameters are given by a^^hl^ where a^^h = I'^e.h^? is 
the Bohr radius for the electron and hole respectively. Taking B ~ 8T and = 12.53 yields 
(GaAs) ahjl = 0.2 and a^/l = 1.1, thereby rendering any such perturbation approach unre- 
liable. For this reason we diagonalize the charged exciton Hamiltonian exactly numerically. 
We perform the diagonalizations in regions of the Hilbert space charaterized by quantum 
numbers for the charged exciton; i.e. M, S and S^. The actual basis used has 5 hole Landau 
levels, 2 electron Landau levels and 16 angular momentum states per Landau level. A larger 
number of hole Landau levels is used because of the small hole cyclotron energy. A finite 
number of angular momentum states is used since the confinement breaks the degeneracy 



of the Landau level. Clearly as B — > an increasingly large number of Landau levels should 
be used. We therefore only quote our calculations for B> 2 T, the reliability of the answers 
increasing rapidly with B-field. From comparison with published experimental data given 



later in this paper and in Ref. [|I^ the authors believe this basis to be sufficient for captur- 
ing the essential physics of the charged exciton. For the purpose of comparison, the exciton 
states are calculated using the same basis in a manner identical to that used for the charged 
exciton. The reliability of the calculation has been verified by reproducing known results 
as follows: (a) in the limit of zero confinement the hidden symmetry result is recovered for 



B— >■ oo, (b) the B=0 confined 2-electron results of Ref. |^ are reproduced, (c) the main 
results found in a previous study on an ideal 2D charged exciton ITH] are obtained. 



III. PHOTOLUMINESCENCE SPECTRA 

In this section we consider the photoluminescence (PL) spectra of charged and uncharged 
excitons. These spectra arise from the recombination of a valence hole with a conduction 
electron giving out a photon of energy huj such that huj = Ei — E j where Ei and Ej are the 
initial and final energies respectively. The recombination process conserves total angular 
momentum. The emitted photon has angular momentum = ±1 which is taken up by 
the atomic part of the Bloch wavefunction (i.e. transition between s ans p orbitals) causing 
the angular momentum of the envelope function to remain unchanged under recombination. 
Throughout this paper we are only dealing with envelope wavefunctions (i.e. we neglect 
bulk K.p bandmixing effects) hence we have the selection rule AM = 0. Consequently, 
only excitons with M=0 can undergo recombination and contribute to the PL spectrum 
(X(M = 0) —* '-/). Now consider recombination from charged excitons. Angular momentum 
must still be conserved, but on recombination an electron is left in the final state. This 
electron can have any angular momentum, thereby allowing charged excitons with total 
angular momentum, Mt^ also to contribute to the PL spectrum via X~{M) — > 7 + e~(M). 
The final electron can be in any Landau level provided it has angular momentum M. The 
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higher the energy of the final state electron, the lower the PL photon energy. In this section 
we always consider electrons to be in the lowest Landau level in the final state. Electrons 
left in excited states (i.e. higher Landau levels) are considered later in this paper when we 
deal with shakeup processes. 

We consider the recombination to involve the ±| heavy holes; this is consistent with 
the relevant experimental work These holes combine with ±| electrons giving two 
circularly polarised lines, cr"*" and a", which are Zeeman split by the magnetic field. The 
spin of the remaining electron from the charged exciton has no effect on this splitting as its 
Zeeman energy makes the same contribution to the energies of the initial and final states. 
Experimentally the Zeeman splitting of the exciton PL lines is found to be less than that 
of the charged exciton lines 0]. For clarity therefore the calculated PL lines presented here 
are situated at the average energies of the and o~ lines (i.e. taking g^.^^ — > 0). 

As well as the PL energies, we are interested in the PL strengths. Ignoring complications 
due to population dynamics, we will take the measure of the PL strengths to be | ^0|I/|X^ p 
and I {e^\L\X^'^ p for X and X~ respectively where L = J ■?/'e(r)'?/5/j(r)c/^r. This leads to 
expressions for PL strengths as follows: 



noMo\L\S = 1,S, = ±1 



noMo|L|^ = 1,^^ = 



noMo\L\S = 0,S, = 



0\L\X 
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where 



pT/S _ V- (^T/S 
^ ~ Z-^ ^noMo,n-M,nM 

noA'Io>n-M 



/ + X! Cn-M,noMo,nM 
noMo<n—M 



■"n-M'nM 



n,M 



Here S is the total electron spin, Sz is the z-component of total electron spin and C^''^ are 
the expansion coefficients of the singlet/triplet charged exciton. The are the expansion 



coefficients of tfie exciton. Tfie electron remaining from the charged exciton recombination 
is in the state with quantum numbers {no, Mq}. It is difficult to predict which of the triplet 
(S=l) states contributes most to the PL spectrum. Therefore we use as a measure of 

the PL strength of the triplet line. 

We now turn to the photoluminescence (PL) spectra of charged and neutral excitons in 
a weak in-plane confinement. At all B-fields examined (i.e. B=2 — > 12 Tesla) we find that 
the ground state of the neutral exciton (X) has total angular momentum M — 0, the ground 
state of the singlet charged exciton {X~) has M — and the ground state of the triplet 
charged exciton {Xf) has M = — 1. These all contribute to the PL spectrum. Other charged 
exciton states with M< can also luminesce via X~{—M) — > 7 + e~(— M); the electron in 
the final state is left in the lowest Landau level. Since these states are very close in energy 
to the absolute X~ and Xf ground states and their PL lines are also similar, they must be 
included in a calculation of the PL spectrum. Increasing \M\ increases the energy of the X~ 
state, so we only need to consider those states for which —3 <M< in our calculations. 

In a weakly confined system it makes sense to discuss whether the X~ species are stable 
to the decomposition X~ ^ X + e~ . A charged exciton will decompose if its initial energy 
is greater than that of the final state consisting of an exciton in its lowest energy state 
and a free electron in the lowest Landau level. In a truly unconfined system, the angular 
momentum of the final state electron is irrelevant due to the degeneracy of the Landau 
level. In our calculations however, the electron single particle energies do vary weakly with 
Mc thus we assume angular momentum conservation in order to take account of this; i.e. 
X~{M) X{0) + e~{M) where the ^(0) state is the ground state of the magnetoexci- 
ton. This decomposition is closely related to the PL recombination process for X~. Since 
X(0) can recombine directly {X — > 7) it is clear that an X~ state will be stable to this 
decomposition provided its PL energy is lower than that from the bare exciton. 

Figure 1 shows the development of the spectrum with B-field for a rod length L — 31 
(240A at B=10T) and a weak in-plane confinement. To facilitate computation, the rod 
length is allowed to scale with the characteristic length 1. Recent work on the interface 
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exciton indeed suggests that some reduction of L with B is expected [|22|. The spectra are 



robust to changes in L and hence Fig.l retains the experimental features simulated in our 
earlier paper |T^ using fixed L = 220A. Reference only fully considered PL spectra 



from X~ species with M = 0. In contrast Fig.l includes optical recombination from the 
lowest energy singlet and triplet states corresponding to M = —3, —2, —1 and 0. The X~ , 
and exciton (X) PL spectra are shown separately. A Gaussian broadening of O.SmeV 
FWHM has been introduced. 

The spectra in Figure 1 agree well with those observed experimentally in quantum wells 
of width 200 - 300 A . The PL energies increase with B in a similar way to the experimental 
data in addition the theoretical spectrum consists of a predominately X~ peak about 
ImeV below the exciton line at 8T (as compared to the experimentally observed value of 
1.7meV at 8T). At low B Fig.l shows an X^ peak above the exciton line. This is due 
to recombination from the M = triplet charged exciton. This will most likely not be 
observed experimentally since its position indicates that it is unstable to decompostion 
into an exciton and electron (as discussed earlier) and as such is not expected to form. 
However, on increasing B, this X^ line can be seen to move through the exciton line, thus 
becoming stable to decomposition and allowing it to be observed. This occurs at around 
6 Tesla. The appearance of an X^ line near the exciton (X) line at finite B is also a 
feature of the experimental spectra in which it occurs at around 2.5T [^]. The lower energy 
contribution to the triplet spectra in Fig.l originates from the X^ states with M< 0. As 
B is increased these clearly lose strength to the M=0 triplet line. This is a precursor to 
the hidden symmetry mentioned earlier; as B increases the low energy states are forced to 
have more weight in the lowest electron and hole Landau levels. Apart from the M< 
triplet lines, the PL strength of the lines generally increases with B. This is much more 
pronounced for the exciton (X) and X^ than for X^. This effect is observed experimentally 
and agrees quite closely with Fig.l. The splitting between the exciton (X) and X~ 
lines in Fig.l increases with B in a manner similar to that observed experimentally and 
is more pronounced than that found in our previous study which used a fixed rod length 
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T7| . This suggests that the observed increased binding of the X~ with respect to the 
exciton is to a certain extent due to the scahng with magnetic length of the extent of the 
exciton and charged exciton wavefunctions perpendicular to the well. We emphasize that the 
spectra presented in Fig.l contain recombination from charged excitons with initial angular 
momenta, M = —3, —2, —1, and 0. The calculated spectra clearly have all the qualitative, 
and many of the quantitative features of the experimental spectra. Thus charged excitons 
with a range of initial angular momenta could be involved in the experimental systems 
reported to date. 

We now consider the PL signatures of charged and uncharged excitons in an intermediate 
in-plane confinement potential (cUg ~ Uc). The ground states are found to have angular 
momentum M = for the singlet charged exciton, M = — 1 for the triplet and M = 
for the exciton. This is identical to the weakly confined case. Again we examine the PL 
spectrum from charged excitons with M< as these undergo recombination leaving the 
remaining electron in a state of M< which becomes part of the lowest Landau level at 
high enough B-field. 

Fig. 2 shows the calculated PL spectra for the lowest energy singlet and triplet X~ states 
with —3 <M< together with the PL line from the ground state of the exciton (M=0). The 
exciton is included for comparison but would not be observable in the same quantum dot 
as a charged exciton since the parabolic confinement holds the electrons and holes together 
thus preventing complete ionization. The same rod length is used as in Fig.l. The in- 
plane spatial extent of the charged excitons ranges from 300A to 400A at B=10T (discussed 
later) which is larger than the rod length, L (220 A at B=10T) ensuring that the charged 
excitons retain the quasi-2D nature required by our model. These spectra contain interesting 
structure. The singlet {X~) spectrum is always below the exciton (X) line and the triplet 
{X^) spectrum has two parts, one higher than the exciton line which corresponds to the 
M=0 triplet line, and one at similar energies to the singlet spectrum which corresponds to 
the M 7^ triplet states. As B is increased the X-X~ splitting increases slightly (1.2meV 
at 2T to 1.5meV at 12 T) in a similar manner to that seen in the weakly confined case. 
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We also note that the major group of hnes show a diamagnetic shift upwards in energy 
with increasing B. More interesting is the development of the singlet and triplet spectra 
with B. We shall first consider the singlet spectra. At low B the only strong line arises 
from the singlet ground state (M = 0). The M = —1, —2, —3 lines are arranged below this 
decreasing in both PL intensity and energy as M decreases. As B is increased these lower 
M lines increase their PL strengths and energies moving through the M = singlet line 
around B=7T. At high B the lower M singlet lines are clearly arranged at higher PL energies 
than the M = line with PL strengths sufficiently high to allow them to be observed. The 
ordering with increasing PL energy is M = 0,-1, —2, —3 at high B. The triplet spectrum 
below the exciton line shows the same development with B- field. As B is increased the 
M = —2, —3 triplet lines gain PL strength and move through the ground state (M = — 1) 
triplet line around B=7T arranging themselves in the order M = — 1, —2, —3 with increasing 
PL energy at high B . The interesting development of this structure with B is due to the 
competition between the effects of the B-field and the parabolic confinement mentioned in 
Sec. y. This competition is very important in this case of intermediate confinement. At low 
B the parabolic in-plane confinement dominates; this removes the degeneracy of the Landau 
levels and vastly lowers the number of states available to the charged exciton with which to 
minimize its electrostatic energy. In the lowest energy classical configurations the hole has 
Mh = while the electrons have M^^ ~ such that Mg^ + M^^ = M, where M is the 
angular momentum of the charged exciton. This configuration has negligible overlap between 
electron and hole single particle wavefunctions giving the very low PL strength at low B for 
higher \M\ as observed in Fig. 2. As B is increased the effect of the confinement begins to 
disappear as the Landau levels form. The degeneracy at a given Me and Mh now allows the 
electrostatic energy to be minimized with finite single particle electron-hole overlap giving 
the stronger PL strength. This competition between B-field and confinement also explains 
the changes in the PL energy orderings. When the confinement dominates (i.e. low B-field) 
the charged excitons with higher \M\ have to leave their extra electron in a higher \M\ (i.e. 
higher energy) state after recombination thus pushing the PL line to lower energies. This 
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effect disappears as B increases and begins to dominate since tfie degeneracy of the Landau 
levels makes the energies of the final M states increasingly similar. It makes sense, therefore, 
that the cross-over between these effects occurs around 6.5T, which is precisely where the 
expansion parameter, 2u!e,h/^^e,h — 1- 

Finally we turn to the case of a strong in-plane confinement. Fig. 3 shows the PL spectra 
for a rod length L = 31 and the strong confinement potential. Again the lowest energy state 
of the singlet charged exciton has total angular momentum M — 0, the triplet M — —1, 
and the exciton M = 0. We consider here photoluminescence from the lowest energy singlet 
and triplet states with angular momenta —3 < M < 0. The exciton PL from the ground 
state (M = 0) is also shown for comparison. However, as with the case of intermediate 
confinement, decomposition (i.e. ionization) of the charged excitons is not possible here 
since the strong confinement keeps the number of particles in the quantum dot constant. 
We note that the spectral hues posess a diamagnetic shift upwards as B increases. The 
form of the spectrum is essentially unchanged over the range of B-fields shown. This arises 
because the confinement dominates at all B-fields shown. Thus we see a spectrum which is 
very similar to the low B- field spectra in Fig. 2. The only hues with significant PL strength 
are the ground state singlet (M = 0), the ground state triplet (M = —1) together with the 
lowest M = triplet state, and the M = exciton ground state. The PL strengths of the 
other lines are low due to the strong confinement as explained above for the intermediate 
confinement case. The M — triplet line lies far above the exciton line. The M — singlet, 
M = —1 triplet and M = exciton lines are grouped within 4meV of each other for all 
B-fields in Fig. 3. The triplet M — —1 line is now the lowest PL line, then comes the singlet 
M — followed by the exciton. It is only with this very strong confinement that the triplet 
line drops below the singlet in the PL spectra. 
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IV. SHAKEUP PROCESSES 



In the previous section we considered recombination from charged excitons in which the 
electron in the final state was left in the lowest Landau level, i.e. rig = and Me < 0. We 
now consider recombinations in which the final electron is left in higher Landau levels. These 
are known as shakeup processes. We will concentrate on cases where the final electron is left 
in the first excited Landau level. These we denote as SUi lines in the PL spectra. There 
are two ways in which an SUi PL line may occur. In both cases the process is X~{M) — »• 
7 + e~(ne, Me = M); one case has = 1, Mg < 0, the other involves rie = 0, Me = 1. We 
label the two processes SUi(n = 1) and SUi(n = 0) respectively. These shakeup processes 
have significantly lower PL strengths than the recombinations dealt with in the previous 
section. They still obey the AM = selection rule, however the dipole matrix element is 
very small. In all Figures showing shakeup processes in this paper the strong PL spectrum 
is also shown for comparison. 

Figures 4(a) and (b) show the SUi(n = 1) and SUi(n = 0) PL lines respectively for the 
case of a weak in-plane confinement with rod length L = 31. These spectra are very similar 



to recent experimental data for SU processes obtained in quantum wells |jT2|. Consider 
Figure 4(a) which shows the SUi(n = 1) processes. The energies agree very well with the 
observed values (eg the splitting between the shakeup X' lines and the non-shakeup X~ 
lines is 7.4meV at 4T in Fig.4(a) and 6meV experimentally |]l2l). As B increases the SU 
lines drop almost linearly. Given the weak confinement we expect the difference between 
the non-shakeup X~ line and the shakeup line to be huj^, the cyclotron frequency associated 
with the remaining excited electron. This is indeed what we find; the SUi(n = 1) line drops 
by 1.3meV/T while the nonshakeup X~ rises by about 0.4meV/T, the difference between 
the two lines being 1.7meV/T corresponding to hu^ as expected. Experimentally the SUi 
line is found to vary linearly with B, however the difference between it and the X~ is 
slightly less than hu^ which has been tentatively attributed to an increase in the electron 
mass due to penetration of the AlGaAs layer . Looking in more detail at the calculated 
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spectrum we see that the singlet SUi(n = 1) hnes are very weak in comparison to the triplet 
SUi(n = 1) lines. The latter dominate the SUi(n = 1) spectrum at low B-field. At B=2T 
the SU spectrum is about 100 times weaker than the strong non-shakeup X~ line, this 
compares very well with the observed intensity ratio from Ref. |T^. As B increases the main 
triplet SU line drops rapidly in PL strength thus rendering the SUi(n = 1) unobservable at 
higher B. This is again found experimentally. The main triplet line consists predominately 
of M = —1, —3 lines, the M = 0, —2 lines being similar in strength to the singlet lines. The 
large drop in shakeup PL strength with B is due to the increased energy difference between 
Landau levels causing the initial charged exciton states to have less contribution from both 
electron and hole higher Landau levels. 

Fig. 4(b) shows the SUi(n = 0) processes. Again the energies agree well with observed 
values |T^. The SU spectrum is in fact almost identical to that of Fig.4(a) except that 



now the singlet line dominates and the triplet line is very weak. Again at B=2T the SU 
line is about 100 times smaller than the strong non-shakeup X~ line. The gradient of 
the SUi(n = 0) line is 1.3meV/T, the same as that in Fig. 4(a). The SUi(n = 0) line loses 
strength rapidly with increasing B-field. This is again due to the increased energy difference 
between Landau levels causing the charged exciton to have less contribution from both 
electron and hole higher Landau levels. While both SUi(n = 1) and SUi(n = 0) processes 
match very well with the observed SUi line, it is likely that the majority of the observed 
PL strength actually comes from the SUi(n = 1) process which is dominated by the M = 
— 1 triplet line. This is the lowest energy triplet state and has a weak PL strength for 
a recombination process which leaves the final electron in the lowest Landau level thus 
suggesting that there will be a large population to undergo recombination via a shakeup 
process. 

Figures 5(a) and (b) show the SUi(n = 1) and SUi(n = 0) lines respectively for the case 
of an intermediate in-plane confinement with rod length L = 31. Again the non-shakeup 
lines are included for comparison. In both Fig. 5 (a) and (b) the PL strengths of the SUi lines 
are significantly lower than the non-shakeup spectrum. First we consider the SUi(n = 1) 
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lines in Fig. 5 (a). The overall SUi spectrum is weaker than for the case of weak confinement 
(cf. Fig.4(a)). The triplet SUi lines dominate over the singlet lines and there is no longer 
a linear drop of PL energy with increasing B-field. Also the splitting between the shakeup 
X~ lines and the non-shakeup X~ lines does not change linearly with B. This is due to the 
confinement which prevents hZJeA from varying linearly with B. The competition between 



confinement and B-field was important in the non-shakeup processes of Sec. |lTl| , spreading 
the PL energies of the lines and affecting their strengths. Here we also see an effect on the 
PL energies of the SUi(n = 1) lines. At low B-field, where the confinement dominates, the 
triplet SU PL lines are separated, the M = —3 state having the lowest PL energy then 
M = —2, —1 and finally M = with the highest energy. The energy separation is due to 
the confinement which raises the energy of the —M final electron state thus lowering the 
PL energy. This effect is removed as Landau levels form with higher B-field (B>6.5T). The 
M=0 triplet SUi line always appears at higher energy because its initial state is higher in 
energy. 

We now consider the SUi(n = 0) processes of Fig.5(b). The singlet spectrum clearly 
dominates the overall shakeup spectrum. It decreases linearly with increasing B and its 
intensity decreases with B also. The energy of the singlet SUi(n = 0) line is higher than 
the SUi(n = 1) lines of Fig.5(a). Finally we notice that the high PL strength of the triplet 
SUi(n = 0) line at B=2T is large, being in fact only ten times weaker than the non-shakeup 
lines. This occurs because the dominant effect of the confinement at low B-field is to force 
the ±Me^h states to have similar energies, thereby causing the M=l triplet charged exciton 
to configure itself more like the M< charged excitons and luminesce accordingly. 

Figures 6(a) and (b) show the SUi(n = 1) and SUi(n = 0) spectra respectively for the 
case of strong in-plane confinement, with a rod length L= 3/. Again the non-shakeup spectra 
are shown for comparison. Considering the SUi(n = 1) processes in Fig. 6(a), we see that 
the triplet lines dominate the singlet lines. The triplet spectrum is clearly split into four 
sets each corresponding to a shakeup process from a charged exciton with a given M. The 
M = —3 has the lowest SU energy followed by M = —2, —1 and finally M = 0. This 
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separation is due to the strong confinement (as in the low field regime of Fig. 5 (a)) which 
raises the energy of the — Mg final state electron thus lowering the PL energy. As B increases, 
the —Me states start to converge in energy to form the Landau level; the spacing between the 
lines is consequently reduced. The SUi(n = 1) lines are very weak in this strong confinement 
hmit. This is because the strong confinement forces the charged exciton to have very little 
contribution from the n 7^ states which now have much higher single particle energies. 

The SUi(n = 0) lines of Fig. 6(b) are much higher in PL energy than those of Fig. 6(a). 
The singlet line decreases linearly with B having a gradient of about ImeV/T. The triplet 
line is not linear, and for low fields is close to the non-shakeup spectrum in both energy 
and PL strength. As in the case of intermediate confinement at low B, this is due to the 
similarity in energy of ±Me,/i states, thereby causing the M — 1 triplet charged exciton to 
configure itself in a similar manner to the M = —1 state and luminesce accordingly. We 
notice this effect dies away with increasing B. However, it clearly persists to higher B than 
in the case of intermediate confinement. 

From Figs 4,5 and 6 we see that the SUi(n = 1) spectra become increasingly disimilar 
from the SUi(n = 0) spectra as the in-plane confinement is increased. This occurs because 
the degeneracy of the Landau levels is lost as the confinement is increased. The single particle 
states instead form themselves into harmonic oscillator levels. Thus for all confinements the 
SUi(n = 1) processes leave the final electron in an excited state, whereas in the SUi(n = 0) 
processes, the final state electron is increasingly in the lowest M=l oscillator state. As such 
the SUi(n = 0) process is more like a normal recombination process than a shakeup process. 

V. ELECTRON-HOLE DISTRIBUTION FUNCTIONS 

Our exact diagonalization of the charged exciton Hamiltonian gives us the full expansion 
coefficients for the various charged exciton wavefunctions. We can use these to investigate 
the structure of the charged exciton states. In particular, we are interested in how the 
electrons are correlated with the valence hole. This is conveniently visualized using the 
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electron-hole pair distribution function, g^r) 



g{r) = J P(r + ro, ro)rf'ro (10) 

where 

P(r,ro) = (|E^(re,-r)5(rh-ro)^ (11) 

Thus g{r) can be interpreted as the electron density relative to the hole and P(r, ro) as the 
probability of finding an electron at r and a hole at Tq. P(r, Fq) is written in terms of our 
previous singlet and triplet wavefunction expansion coefficients {C^ and C^) as follows: 

^^/^(r,ro)=Ehr±5rr (12) 

nm 

where S takes the positive sign, T takes the negative sign and A and B are given by: 



nimi>nm 
nm>n2m2 

^nm ~ ^nm ~^ ^nm..nm,.n,.m.,, nhrrih (^o) ^ram (l") 

^nm — ^nm ^ (l^) 

The pair distribution functions are radially symmetric so we are able to plot the radial 
pair distribution function (RPDF) given by p(r) = 27rr5f(|r|). Fig.7 shows how these RPDFs 
are used to determine the structure of the charged exciton. If both electrons are equally 
correlated with the hole (Fig. 7(a)) then we expect p(r) to have a single peak as in Fig.7(b). 
If however one of the electrons is closer to the hole than the other, as in Fig. 7(c), then we 
expect to see two peaks corresponding to the inner and outer electron rings as in Fig. 7(d). 
Charged excitons with this doubly peaked RPDF have a structure which is more like an 
exciton with an orbiting outer electron {X + e~) 

We now examine in more detail the electron-hole radial pair distribution functions for 
a charged exciton with rod length L= 3/. Figure 8 shows how the electron-hole RPDF 
for the lowest energy singlet state (M=0) with a weak in-plane confinement varies with 
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B-field. As B increases the electrons become more closely correlated with the hole. As is 
generally the case, the overall form of the functions doesn't change; instead they just scale 
with the characteristic length /. In the case of weak confinement, this length is dominated 
by the magnetic length giving rise to a strong scaling with B. For stronger confinements the 
characteristic length becomes more dependent on the confining potential, thereby reducing 
the effects of changing B-field on the distribution functions. 

Figures 9 and 10 demonstrate how the electron-hole RPDFs of the lowest energy singlet 
and triplet states for various M depend on the confining potential. Figures 9(a) and (b) 
are for weak in-plane confinement, at B=10T, with singlet and triplet charged excitons 
shown respectively. Figure 9(a) exhibits two different types of curve. The M < curves 
are fairly narrow curves peaked at around 100 A indicating that both electrons are equally 
close to the hole. The M >0 curves are broader, peaking at between 200 and 300A. The 
M — 3 curve even shows a shoulder at higher electron-hole separation indicating an X + e~' 
structure. Figure 9(b) shows the corresponding triplet curves. Again there are two types of 
curve. The triplet M < electron-hole pair distributions are narrow and peaked at about 
150A, whereas the M > curves are much broader. As M increases from to 3 a distinct 
shoulder develops at higher electron-hole separation indicating the X + e" structure. In 
general, the triplet curves broaden more quickly than the singlet curves with increasing M. 
Of particular interest is the M — triplet electron-hole RPDF. This shows a broadening 
that the corresponding singlet curve does not. The reason for this arises from the symmetries 
of the wavefunctions. The singlet charged exciton has electrons which are symmetric under 
interchange of their spatial wavefunctions. Thus the M — singlet charged exciton can have 
both electrons and the hole with the same angular momentum Me,h — 0. This gives rise to 
the sharply peaked singlet M = curve in Fig.9(b). The triplet charged exciton however 
has a spatial wavefunction which is antisymmetric under electron interchange. This implies 
that the lowest energy M — triplet state cannot be built with both electrons and the hole 
in states with M^^h = 0. Instead the state will have one electron with = 0, the other 
with Me = —M and the hole with Mh — M. A state with Mh < M would force the electron 
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with Me = to change to a state of higher Mg, thus placing it in a higher Landau level; this 
is energetically unfavourable. A state with > M would force the hole into a higher \M\ 
state than either electron; the electrons would hence come closer in Mg. This situation is 
clearly electrostatically unfavourable given the relation between \M\ and < > provided 
by Eqn. ^. Thus from symmetry and energy considerations, the lowest energy M = triplet 
charged exciton should have one electron closer to the hole than the other. This gives rise to 
the broad electron-hole RPDF seen in Fig.9(b). It also partially explains the fact that the 
triplet ground state has angular momentum M = — 1, whereas the singlet ground state has 
M = 0. The electron- hole RPDFs at B=10T in the case of strong confinement are shown 
in Fig. 10. As with Fig.9, only the lowest energy states at each M are considered. Figure 
10(a) shows the electron-hole RPDF for the singlet charged exciton. The M = case is still 
the most sharply peaked. The other RPDFs all have similar peak heights and widths. This 
is a direct result of the strong confinement which lifts the degeneracy of the Landau levels 
instead giving instead single particle states corresponding to a harmonic oscillator. Now 
the degeneracy is between single particle states with the same |M|. This forces the —M 
charged exciton states to behave very like those with +M. This can be seen in the RPDFs of 
Fig. 10(a) where the M = — 1 curve resembles the M = +1 curve etc. Comparing Figs. 10(a) 
and 9(a) we see that the distribution functions peak at similar electron- hole separations 
in both cases. The high confinement is conspicuous however in removing the shoulders 
corresponding to high electron-hole separation from the distribution functions; consequently 
all these lowest energy states have both their electrons equally well correlated with the hole. 
Figure 10(b) shows the corresponding electron-hole RPDFs for the triplet charged excitons. 
This shows strong similarities with the singlet case. The main distinction is the sharply 
peaked triplet M = — 1 RPDF which is narrower than the others. This corresponds to the 
M = RPDF in the singlet case (Fig. 10(a)). It comes as no surprise therefore that the 
M = state is the lowest energy singlet state whilst the lowest energy triplet state has 
M = —1. As in the weakly confined case this difference in the angular momenta of the 
lowest energy singlet and triplet states is due to the symmetries of the singlet and triplet 
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spatial wavefunctions under electron interchange. Again we see the triplet states showing 
evidence of underlying oscillator single particle states since the +M and —M states have 
become very similar. 

The electron-hole RPDFs for an intermediate in-plane confinement behave as the strong 
or weak confinement cases depending whether the B- field or the confinement is dominant. 
Earlier we saw that the cross-over between the dominance of the confinement or the B-field 
occurs at about B=6.5T for this intermediate confinement. Below B=6.5T the electron-hole 
RPDFs are similar to those in Fig. 10. Above B=6.5T the electron-hole RPDFs are similar 
to those in Fig. 9. 

We now turn to the excited states. As an illustration we look at the excited states with 
the same M as the lowest energy singlet (M = 0) and triplet (M = — 1) charged excitons. 
The electron-hole RPDFs for the first five states of these systems are shown in Figs. 11 (a) 
and (b). A charged exciton can become excited in essentially two different ways. First it 
can start to ionize itself by moving the second electron increasingly further away from the 
electron and hole which are left in the centre. Second it can keep both electrons equally 
close to the hole, but alter their correlations with the hole. In the first case we should see 
two peaks forming in the electron-hole distribution function. In the second the distribution 
function should have a single peak, but the behaviour near the origin and the peak position 
should be modified . With this in mind we consider Fig. 11. First we consider the singlet case 
(Fig. 11 (a)). The lowest four energy states clearly show the electron- hole RPDF becoming 
broader as the energy increases, until eventually at level 4 there is a definite high electron- 
hole separation shoulder indicating a partial separation into X + e~ . The next energy level 
up sees both electrons equally correlated with the hole, the RPDF vanishing more quickly 
at the origin and its peak position having moved to higher separations. This corresponds 
to the second type of excitation. The triplet energy levels at M = — 1 (Fig. 11(b)) show a 
very similar behaviour. The electron-hole RPDFs broaden as the energy level is increased. 
A higher separation shoulder can be just distinguished in the RPDF for level 5. 

Similar effects are seen in the excited states at different M and different confinements. 
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The clearest case for showing the two types of excitation process is a strongly confined 
charged exciton with M = —2 at B=10T. This is shown in Fig. 12. Here the first three 
levels show broad electron-hole RPDFs. Level 4 has two distinct peaks in its RPDF arising 
from two different electron rings around the hole. Level 5 then has an electron-hole RPDF 
which has one peak moved to slightly higher e-h separation, and clearly has a different 
behaviour at the origin. This suggests that the process of excitation for a charged exciton 
at fixed angular momentum M involves repetition of the two-step process corresponding to 
Figs. 7(b) and 7(d). Each cycle pushes the peaks to higher r and forces the RPDF to vanish 
more quickly at the origin. 



VI. SUMMARY 

To summarize, we have studied a negatively charged exciton in a uniform B-field by 
means of numerical diagonalization. We have discussed our results in terms of the PL spectra 
for various in-plane confinements and experimentally relevant B-fields. The PL spectra were 
calculated considering recombination from a range of initial angular momentum states of the 
charged exciton (—3 < M < 0). A realistic form factor consistent with 200-300A quantum 
wells has been used to modify the interactions. We found that our weak in-plane confinement 
results agree very closely with published experimental spectra for high mobility quantum 
wells (Ref. p). We also get good agreement with the PL energies and intensities of published 
shakeup lines in the PL spectra (Ref. |]12|)- These results demonstrate that the expected PL 
signatures of charged excitons are not qualitatively affected by the presence of a range of 
initial angular momentum states. The one exception to this is the triplet line which appears 
at finite B-field. This is due to the M = triplet charged exciton which must be present for 
the line to be observed. We also discussed the PL spectra (including shakeup lines) for the 
cases of intermediate and strong in-plane confinements. These represent predictions of the 
expected PL spectra from charged excitons in quantum dots. Our calculated spectra show 
interesting structure which should be observable and we hope this will stimulate further 
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experimental study on such systems. Finally, we calculated various electron-hole radial pair 
distribution functions for the charged excitons. These allowed us to examine in detail the 
structure of the charged exciton states as a function of angular momentum, B-field and in- 
plane confinement. The excited states of the charged exciton were also studied. We found 
that excitation occurs via two methods. First the charged exciton can undergo partial 
ionization whereby one of the electrons moves away from the hole leaving a configuration of 
an exciton orbited by an outer electron. Second, both electrons remain equally close to the 
hole, but their correlations with it change. The combination of these two processes provides 
insight into the excitation mechanism for charged excitons at a given angular momentum. 

We would like to thank A.J. Shields and G.Finkelstein for preprints. Funding was pro- 
vided by EPSRC through direct support (J.R.C.) and through Research Grant No. GR/K 
15619 (N.F.J, and V.N.N.). 
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Figure Captions 

Figure 1: PL Spectra for a weak in-plane confinement and rod lengtfi L — 31 (240A at 
B=10T). Recombination from the lowest energy singlet {X~) and triplet (Xf) charged 
excitons with M = —3, —2, —1, is included. Recombination from the exciton (X) ground 
state (M = 0) is also shown. A Gaussian broadening of 0.3meV FWHM is included. 

Figure 2: PL Spectra for an intermediate in-plane confinement and rod length L — 31. 
Recombination from the lowest energy singlet {X~) and triplet (X^) charged excitons with 
M = —3, —2, —1, is included. Recombination from the exciton (X) ground state (M = 0) 
is shown for reference. A Gaussian broadening of O.SmeV FWHM is included. 

Figure 3: PL Spectra for a strong in-plane confinement and rod length L — 31. Re- 
combination from the lowest energy singlet {X~) and triplet (X^) charged excitons with 
M = —3, —2, —1, is included. Recombination from the exciton {X) ground state (M = 0) 
is shown for reference. A Gaussian broadening of 0.3meV FWHM is included. 

Figure 4: Shakeup lines in the PL spectrum for a weak in-plane confinement and rod length 
L — 31. (a) SUi(n = 1) shakeup from the lowest energy M — —3, —2, —1,0 singlet and 
triplet charged exciton states, (b) SUi(n = 0) shakeup from lowest energy M = 1 singlet 
and triplet charged exciton states. The non-shakeup spectra of Fig.l are also shown for 
reference. A Gaussian broadening of 0.3meV FWHM is included. 

Figure 5: Shakeup lines in the PL spectrum for an intermediate in-plane confinement and 
rod length L = 31. (a) SUi(n = 1) shakeup from the lowest energy M = —3, —2, —1, singlet 
and triplet charged exciton states, (b) SUi(n = 0) shakeup from lowest energy M — 1 singlet 
and triplet charged exciton states; note that different magnifications are used for the singlet 
and triplet processes. The non-shakcup spectra of Fig.2 are also shown for reference. A 
Gaussian broadening of 0.3meV FWHM is included. 

Figure 6: Shakeup lines in the PL spectrum for a strong in- plane confinement and rod 
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length L — 31. (a) SUi(n = 1) shakeup from the lowest energy M — —3, —2, —1, singlet 
and triplet charged exciton states, (b) SUi(n = 0) shakeup from lowest energy M = 1 singlet 
and triplet charged exciton states; note that different magnifications are used for the singlet 
and triplet processes. The non-shakeup spectra of Fig.3 are also shown for reference. A 
Gaussian broadening of O.SmeV FWHM is included. 

Figure 7: Schematic diagram showing radial pair distribution functions (RPDF) for typical 
configurations (a) Both electrons equally correlated with the hole, (b) The electron-hole 
radial pair distriution function p(r) associated with part (a), (c) One electron closer to the 
hole than the other, (d) The electron-hole radial pair distribution function p(r) associated 
with part (c). 

Figure 8: Electron-hole radial pair distribution function p(r) at various B-field for lowest 
energy M — singlet charged exciton in a weak in- plane confinement and rod length L — 31. 

Figure 9: Electron- hole radial pair distribution function p(r) for lowest energy charged 
exciton states with B=10T in a weak in-plane confinement and rod length L = 31. Charged 
exciton states with M = — 3 — > 3 are considered, a) Singlet states, b) Triplet states. 

Figure 10: Electron-hole radial pair distribution function p(r) for lowest energy charged 
exciton states with B=10T in a strong in-plane confinement and rod length L — 31. Charged 
exciton states with M = —3 — > 3 are considered, a) Singlet states, b) Triplet states. 

Figure 11: Electron-hole radial pair distribution function p(r) for the lowest five energy 
states of (a,) M — singlet charged exciton (level 1 is the singlet ground state) and (b) 
M — —1 triplet charged exciton (level 1 is the triplet ground state). B=10T with a weak 
in-plane confinement and a rod length L = 31. 

Figure 12: Electron-hole radial pair distribution function p(r) for the lowest five energy 
states of an M = — 2 singlet charged exciton in a strong in- plane confinement with B=10T 
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and a rod length L — 31. 
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